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The quantum deformation of the Poisson bracket is the Moyal bracket. We construct quantum 
deformation of the Dirac bracket for systems which admit global symplectic basis for constraint 
functions. Equivalently, it can be considered as an extension of the Moyal bracket to second-class 
constraints systems and to gauge-invariant systems which become second class when gauge-fixing 
conditions are imposed. 
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I. INTRODUCTION 



The association rules between real functions in the phase space of unconstrained classical systems and Hermitian 
operators in the Hilbert space of the corresponding quantum mechanical systems are discussed since a long time [1-5] . 

0^ ' The commonly used association rule proposed by Weyl [1] consists in replacing products of the canonical variables 
with symmetrized products of operators of the canonical variables. The Wigner function [2] is associated to the 

C ■ density matrix. 

■ To the lowest order in the Planck's constant, the product fg of operators associated to functions / and g in the phase 
space corresponds to the pointwise product fg. A deformation of the pointwise product which keeps the association 
rule to all orders in the Planck's constant is constructed by Groenewold [31. This product is known as the star-product. 
The Moyal bracket [4] appears as a skew-symmetric part of the star-product. It defines the representation of a 
| commutator — i/H[, ] in the space of the functions. The Moyal bracket is antisymmetric, coincides with the Poisson 
'-^ ' bracket to the lowest order in the Planck's constant, satisfies the Jacoby identity, and keeps the association rule. 

The bracket satisfying these properties is essentially unique [5] . The Moyal bracket governs the quantum evolution of 
systems in the phase space like the Poisson bracket governs the classical evolution. A survey of the star-product and 
the Moyal quantization can be found in [6, 7]. 

_ _ i 

The quantum dynamics of unconstrained systems can therefore be formulated in the phase space in terms of the 
Hamiltonian and Wigner functions, with the pointwise product replaced by the star-product. The average values of 
quantum observables can be computed by averaging the symbols of the Hermitian operators over the Wigner function. 

The specific feature of gauge theories is the occurrence of constraints which restrict the phase space of gauge- 
invariant systems to a submanifold. A systematic Hamiltonian approach to gauge theories and general constraint 
systems and the corresponding operator quantization schemes were developed by Dirac [8] . 

Schemes based on the path integral method have also been proposed and found to be useful for quantization of 
gauge theories [9, 10] and second-class constraints systems [11, 12]. 

Gauge systems are quantized by imposing gauge-fixing conditions which convert them into second-class constraints 
systems. Anomalous gauge theories [13-17], the 0(n) non-linear sigma model [18-20], many-body theories involving 
collective and elementary degrees of freedom [21-23] are second class from the start. 

If constraint equations are solved, the system can be restricted to the constraint submanifold and treated accordingly 
as an unconstrained system. In many cases, however, it is not possible to solve constraint equations. The method 
proposed by Dirac for classical second-class constraints systems solves the evolution problem by applying the Dirac 
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bracket to functions of canonical variables in the unconstrained phase space. It allows thereby to avoid typical 
complications connected to restriction of the systems to constraint submanifolds. 

In this paper, we discuss second-class constraints systems from the Moyal quantization perspective. The main 
problem we focus attention to is the construction of quantum deformation of the Dirac bracket able to govern the 
evolution of quantum constraint systems in the unconstrained phase space. The guiding idea in development of the 
quantization scheme is that many functions in the unconstrained phase space correspond to one and the same physical 
observable. We thus come naturally to the notion of equivalence classes of functions, and also operators. In the next 
Sect., we specify equivalence classes of real functions in the unconstrained phase space and, in Sect. Ill, of Hcrmitian 
operators in the Hilbcrt space. An association rule between equivalence classes of operators and functions, based on 
the Weyl's association rule, is established. In Appendix A, some useful properties are derived for the skew-gradient 
projection in terms of which the equivalence classes of functions and operators are defined. In Sect. IV, we give a 
summary and, in Appendix B, proofs of the basic properties of the Weyl's association rule and of the star-product. In 
Sect. V, the quantum deformation of the Dirac bracket is constructed for systems which admit global symplectic basis 
for constraint functions. In Appendix C, we provide an explicit form of the lowest order 0(h 2 ) quantum correction to 
the Dirac bracket. As an application, in Sect. VI, the quantum deformation of the Dirac bracket is used to formulate, 
in the unconstrained phase space, an evolution equation for the Wigner function of an n — 1-dimensional spherical 
pendulum which represents a mechanical counterpart of the 0{n) non-linear sigma model. 



II. CLASSICAL SECOND-CLASS CONSTRAINTS SYSTEMS IN THE PHASE SPACE 



Second-class constraints Q a = with a = l,...,2m in an 2n-dimensional unconstrained phase space 
((f) 1 , (f> n , 7T 1 , 7r") have the Poisson bracket relations which form a non-degenerate 2m x 2m matrix 



dct{g a ,g b }^o. 



(11.1) 



Two sets of the constraint functions are equivalent if they describe the same constraint submanifold. One can make 
therefore non-degenerate transformations on the constraint functions without changing the dynamics. 

For an arbitrary given point of the constraint submanifold, there is a neighborhood where one may find the equivalent 
constraint functions in terms of which the Poisson bracket relations look like 



{Go,, Gb} — lab 



where 



lab = 



E n 
-E m 



(II.2) 



(II.3) 



with E m being the unity m x m matrix, l a b^bc = —Sac- The upper and lower indices of vectors are discriminated 
according to the rules % = I a bT b , T a = I ab %, l ab = -l ab , {G a ,Gb} = <5£, etc. The scalar product T a y a is invariant 
with respect to the group of linear symplectic transformations Sp(2m). 

The global symplectic basis (II. 2) for constraint functions exists obviously for m = 1 and, also, for systems of 
point particles under second-class holonomic constraints [20] and second-class non-holonomic constraints satisfying 
the Frobenius' condition [24]. The global existence of the basis (II. 2) is proved for systems with one primary constraint 
[25] and for a broader set of systems using additional assumptions [25, 26]. In general case, the global existence of 
the symplectic basis for constraint functions is an opened question. 

The basis (II. 2) always exists locally, i.e., in a finite ncghborhood of any point of the constraint submanifold [24-26]. 
This is sufficient for needs of the perturbation theory. The formalism presented in this work can therefore to be used 
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to fomulate, in the sense of the perturbation theory, the evolution problem of any second-class constraints system in 
the unconstrained phase space. Eqs.(II.2) follow from the local existence of the standard canonical coordinate system 
[27] where Q a with a — 1, m play the role of the first canonical coordinates and Q a with a = m + 1, 2m play the 
role of the first canonical momenta. 

Let us construct skew-gradient projections £ s (£) of the canonical variables £ onto the constraint submanifold £? a (0 = 
using phase flows generated by the constraint functions. From equations 

{MQ,G a (0} = o (n.4) 

using the symplectic basis (II. 2) for the constraints and expanding 

UO = i + x a g a + \x ab g a g b + ... (ii.5) 

in the power series of Q a , one gets 

UO = ^^{-A{^S ai },S a '},...S a '-}g ai ga a ...Qa k . (n.6) 

One can show (see Appendix A) that any function /(£) projected onto the constraint submanifold 

oo 1 

fs(O^T,M^f(O,Q ai },Q a2 },---0 ak }0a 1 g a2 ...ga k (n.7) 

k=0 

satisfies 

/.(0 = /(&(0)- (H.8) 

Eq.(II.4) tells us that variations of £ s (£) along the phase flows generated by the constraint functions G a (Q are zero. 
It means that £ s (£) belongs to the constraint submanifold 

Q a (UO) = 0- (n.9) 

For any function /(£), one gets {£/ a (£), /(£s(£))} = 0- The reciprocal statement is also true: The coordinates on the 
constraint submanifold can be parameterized by £ s . The coordinates describing shifts from the constraint submanifold 
can be parameterized by g a . The functions / can be presented by / — f(£ s , g a ). If / is identically in involution with 
g a , it depends on £ s only. This can be summarized by 

{Ga,f} = 0"f = f(UO)- (H.10) 

Eqs.(II.9) and (11.10) become selfevident if one works in the standard canonical coordinate system. 
An average of a function /(£) is calculated using the probability density distribution p(£) and the Liouville measure 
restricted to the constraint submanifold [28]: 

/Jin t 2m 
J -j-(2nrl[5(g a (0)f(0p(0- (11.11) 

On the constraint submanifold £ s (£) = so f(^) and p(^) can be replaced with / s (0 and p s (Q- 
There exist therefore equivalence classes of functions [29] in the unconstrained phase space: 

f (0 ~ 9(0 " fs(Q = 9.(0- (H.12) 

Here, /(^) ~ g(£) means that the functions are equal in the weak sense, /(£) w g(^), i.e., on the constraint submanifold. 
We shall see that the symbols ~ and « acquire distinct meaning upon quantization. Note that f(£) ~ / s (0- Eqs.(II.8) 
and (II. 9) imply g a ~ 0. 
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Given the hamiltonian function H, the evolution of a function / is described using the Dirac bracket 

^f = {f,H} D . (11.13) 
In the symplectic basis (II. 2), the Dirac bracket looks like 

{f,g} D = {f,g} + {f,G a }{9a,g}- (n.14) 

On the constraint submanifold, one has 

{f,g}n = {f,g s } = {f s ,g} = {f s ,g s }. (n.15) 

Two hamiltonian functions are equivalent if they generate within the constraint submanifold identical phase flows. 
The components of the hamiltonian phase flow, which belong to a subspace spanned at the constraint submanifold 
by phase flows of the constraint functions, do not affect the dynamics and could be different. H and H s are thereby 
equivalent, so Eq.(II.12) characterizes an equivalence class for the hamiltonian functions either. Among functions of 
this class, TL S is the one whose phase flow is skew-orthogonal to phase flows of the constraint functions. 

Replacing H with H s , one can rewrite the evolution equation in terms of the Poisson bracket (cf. Eq.(II.13)): 

p = {f,n s }. (n.16) 

Eq.(II.16) applied to f(t) and g(t) which belong at t = to the same equivalence class, provides with the use 
of Eq.(A.3) and the initial condition / s (0) = g s (0), f s (t) = g s {t). The equivalence relations (H-12) are therefore 
preserved during the evolution. 

If we denote the equivalence class of a function /(£) as £f, the sum of two equivalence classes £f and £ g can be 
defined as £f + £ g = £f+ g on line with f s + g s = (f + g) s , the associative product can be identified with £f g , while 
the skew-symmetric Dirac bracket can be defined as {£f,£ g }o = £{f,g} D - These operations satisfy the Leibniz' law, 

{£f£ g , £h}o = £f{£gi £h}o + {£f, £h}D£ g , (H-17) 

and the Jacoby identity, 

{{£ f , £ g } D , £ h } D + {{£ g , £ h } D , £ f , } D + {{£ h , £ f } Dl £ g } D = 0. (11.18) 

The associative product does not depend, in virtue of Eq.(A.2), on the choice of representatives of the equivalence 
classes. The Dirac bracket can be calculated for arbitrary representatives of the equivalence classes either. Indeed, 
/ ~ g implies 

{/, h} D ~ ({/, h} D ) s = {f s , h s } D = {g s , h s } D = {{g, h} D ) s - {g, h} D , 

where the use is made of Eq.(A.4). 

The physical observables in second-class constraints systems are associated with the equivalence classes of real 
functions in the unconstrained phase space. The equivalence classes £f constitute a vector space O equipped with two 
multiplication operations, the associative pointwise product and the skew-symmetric Dirac bracket {, }d, which confer 
O a Poisson algebra structure. 

The one-to-one mapping £f f s induces a Poisson algebra structure on the vector space of projected functions. 
The sum £f + £ g converts to f s + g s , the associative product £f£ g converts to the pointwise product f s g s , while the 
Dirac bracket {£f,£ g }D becomes the Poisson bracket (cf. transition from (11.13) to (11.16)): 

{f s ,gs}D = {fs,g s }. (11.19) 

These operations satisfy the Leibniz' law and the Jacoby identity and, since (f s + g s ) s = f s + gs, (/sfl's)s — fsgs, and 
{f s , h s } s = {fs, h s } (cf. Eqs.(A.2) and (A. 3)), keep the vector space of projected functions closed. 
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III. QUANTUM SECOND-CLASS CONSTRAINTS SYSTEMS IN THE HILBERT SPACE 

The systems are quantized by the algebra mapping — > f and {,} — > —i/h[,]. To any function /(£) in the 
unconstrained phase space one may associate an operator f in the corresponding Hilbert space. In particular the 
hamiltonian function H(£) and the constraint functions Ga(Q correspond to the operators f) and <8 a , respectively. By 
this mapping the quantal image is called operator. The reverse mapping associates to an operator f a symbol /(£). 

Eqs.(II.2) become 

[<& a ,<d b } = ifiX ab . (mi) 

The Weyl's association rule applied to Eqs.(II.2) yields Eqs.(III.l), provided the quantization is performed in the 
standard canonical coordinate system. 

Any operator in the Hilbert space can be represented as a function of 2n operators f associated to 2n canonical 
variables . 

Let us construct 2n operators f s associated to the projected variables (II. 6). They commute with the constraints 

[Ts,&a]=0. (HI.2) 

The analogue of Eq.(II.6) looks like 

fe=0 

Among the operators acting in the Hilbert space, one expects that, equivalence classes exist. For an arbitrary 
operator f, a projected operator f s can be constructed as follows: 

f» = E T^-[[ f < <8 ° 1 ]. e 02 ]- -& ak ]&a 1 &a 2 ...&a k . (IH.4) 

fc=0 

One has 

[U,e a ] = o (IH.5) 

and 

(f fls ) s = (j.a). = fs9s- (IH.6) 

Two operators f and 9 belong to the same equivalence class provided f s = Q s , i.e., 

f~0^f s =fl s . (III.7) 

The Dirac's quantization method of second-class constraints systems [8] consists in constructing operators repro- 
ducing the Dirac bracket for canonical variables and taking constraints to be operator equations. In the classical 
limit, the commutators of operators (III. 3) satisfy the Dirac bracket relations for canonical variables on the constraint 
submanifold. Furthermore, (<25 a ) s = 0, and so (3 a ~ 0. 

As a consequence of the Jacoby identity and Eq.(III.l), the operator 

= [...[[ ,0 Ol ],0 O2 ],...0°*] (III.8) 

entering Eqs.(III.3) and (III.4) is symmetric with respect to permutations of ai, a^. Any contraction of the upper 
indices with T a b annihilates (III. 8). It follows that 

[C ai - ai - ak ,e ai ] = 0. (III.9) 
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The position of the constraint operators with the lower indices and of the operator (III. 8) in Eqs.(III.3) and (III. 4) is 
not important. One can place, e.g., & ai on the first position, on the second position, etc. 

In order to calculate the average value of an operator, one has to construct the quantal image of the delta functions 
product entering Eq.(II.ll). The projection operator can be written as follows 

Wll-P(^A Q ). (111.10) 

a— 1 

In the classical limit, one recovers the product of the delta functions. 

Let us chose a basis in the Hilbert space in which the first m constraint operators are diagonal, 

e a \g,9. >=g a \g,g* >, (ni.n) 

for a = l,...,m. Congruous to this equation, <8 a might be taken as momentum operators. The last m constraint 
operators can be treated as quantal coordinates. The additional n — m eigenvalues are denoted by g*. 
The projection operator *}3 acts as follows: 

V\g,9* >= \0,9* > ■ (HI.12) 
To arrive at this equation, we split © a = (P A , — Q A ) and A a = (A^, A^) and write 



-— -Y[eM^ a K)\g,g*> 

a—1 



I 



r im \ ft " L ■ '"■ 

dmX '(^ II ex P (|P^A' A ) J] ex P (-lQ-A^)|,,^ > 

V ' A=l B=l 



«^nex P (lp-A^)|. 9 -A", 5 ,> 

/ im \ll m ■ 
^A'^p n ex P (-(/ - >l' A )>l A )\g - X",g* > 

/m 
d m \"l[6(g A -\ A )\g-\",g*> 



A=l 

= \0,g* > 
The average value of an operator f, 

< f >= Trm s v s ] (111.13) 
where r = r + is the density matrix, can be transformed to give 

/rl m nrl n ~ m n 
(Lfi)" <9>9*\W>*e\g,g.> 

° 2wh)n 9 * < S>*lf.t.|0, ff . > = j j^S^ < O lfl .|f.r.|0, fl . >. (IIL14) 

The average values are determined by the physical subspace of the Hilbert space, spanned by the vectors [0, g* >. 
These vectors satisfy the equation 

© a |0,.g*>=0 (111.15) 

which can be recognized as the Dirac's supplementary condition [8] of an equivalent gauge system [20, 24-26], where 
6° with a = 1, m are gauge generators and a with a = m + 1, ...,2m are gauge-fixing operators. 



7 



All density matrices from a given equivalence class correspond to a single physical state, while operators from the 
same equivalence class have equal average values. 

The quantum evolution equation can be written as an extension of the classical evolution equation (11.16) 

ihj t f=[f,Sj s ] (111.16) 

where Sj s is the projection (III. 4) of the Hamiltonian 9). The evolution has the property that at any time f(i) ~ g(£) 
if f(0) ~ g(0). This is suggested by the equation 

[f,fl s ] s = [f s ,fl s ] (111.17) 
which is evident due to (f + q) s = f s + q s and Eq.(III.6). 



IV. WEYL'S ASSOCIATION RULE AND THE STAR-PRODUCT 



The Weyl's association rule can be formulated in a comprehensive way in terms of the operator function [30] 



58(77) = exp(-77 fe j: fe ) 



whose Fourier transform, 



has the properties (see Appendix B) 



/H 2n n i 



<B(0 + 

Tr[®(0] = 1 
d 2n £ 



Tr[»(0®«')] 



f, 

(27rR) n J 2 "^-0. 



Here, 



V - T kl 8 8 

is the so-called Poisson operator. The matrix I kl looks similarly to the matrix (II. 3) 



Tkl 



-E n 
E n 



with E n being the n x n identity matrix. 
Equations 



f(0 = Tr[«8(0fl, 
d 2n £, 



(IV.l) 

(IV.2) 

(IV.3) 
(IV.4) 

(IV.5) 

(IV.6) 
(IV.7) 
(IV.8) 



(IV.9) 



(IV.10) 
(IV.ll) 
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define the Weyl's association rule. 

The canonical variables appear as symbols of operators of the coordinates and momenta: £ 4 = Tr [03 (£)?:*]. 
The phase space of quantum systems is equipped with the Groenewold star-product [3] . Given two functions 

/(0 = Tr[B(Ofl, 
g(0 = Tr[B(Ofl], 

one can construct a third function 

/(O*0(O=Tr[B(Ofg]. (IV.12) 
The star-product is an associative operation. It splits into symmetric and antisymmetric parts: 

f*g = f°g + ^fAg- (iv.i3) 

The explicit form of the star-product can be derived from Eqs.(IV.3) - (IV. 8): 

/(0*ff(0 = /(Oexp(yPMO, ( IV - 14 ) 

where V = and therefore 

f(0°g(0 = /(Ocos(^MO, (iv.i5) 

/(£) A 9(0 = /(0fsin(^) fl (0. (IV.16) 

The Planck's constant 7i appears as a quantum deformation parameter. The antisymmetric product /(£) A is 
known under the name of Moyal bracket. The classical limit of the Moyal bracket is the Poisson bracket: 

hm /(£) A = {/(£), <?(£)}• (IV-17) 

h— >o 

The star-product and the Moyal bracket obey the Leibniz' law 

/ A (g *h) = (/ A g)*h + g* (/ A h). (IV.18) 

This equation is valid separately for symmetric and antisymmetric parts of the star-product. In the last case, 
Eq.(IV.18) provides the Jacoby identity. 

V. QUANTUM SECOND-CLASS CONSTRAINTS SYSTEMS IN THE PHASE SPACE 

The phase space of quantum systems is endowed with the star-product. The corresponding Hamiltonian H(^) and 
constraint functions G a (£) appear as symbols of the operators io and a , respectively, as prescribed by Eq.(IV.lO). 
If the quantization is done in the standard canonical coordinate system, then G a (£) = G a (0- I n general case, the 
equality holds in the classical limit, i.e., 

HmG o (0=aa(0- (V.l) 

Classical systems represent a limiting case of quantum systems, not vice versa, so many quantum systems with the 
same classical limit usually exist. Some ambiguities arise due to the so-called "operator ordering problem". There 
exist, as a consequence, several association rules the most popular of which is the Weyl's rule. The explicit forms of 
distinct association rules are provided by Mehta [5] . The second-class constraints systems have additional ambiguities 
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connected to the choice of constraint functions G a (£). Here, we require only that H(£) and G (£) exist and tend, 
respectively, to and G a (£) as ft ^ 0. 

The analogue of Eqs.(II.2) and (III.l), which specifies the symplectic basis for the constraint functions G ai has the 
form 

G a (0AG b (0=l ab . (V.2) 

These equations are automatically fulfilled provided the constraint operators <3 a obey Eqs. (III.l). 
The skew-gradient projections, analogous to Eqs. (II. 4) and (III. 2), arc defined by 

6(0AG o (0=0. (V.3) 

The projected canonical variables have the form 

oo ^ 

&® = E A Gai) A G<12) - hGak)o Gai ° Ga >- ° Ga - (v - 4) 

fe=0 

(cf. Eqs.(II.6) and (III.3)). The analogue of Eqs.(II.7) and (III.4) is 

OO w 

= E y\ A GQ1) A Ga2) - A ° ak) ° Gai Ga2 - Ga - (v - 5) 

fe=0 

The operation o is not associative in general. The order in which it acts is, however, not important due to commuta- 
tivity of the operators in expressions (III. 3) and (III. 4). In the classical limit one has 

lim /* = /,. (V.6) 

The terms ... o G a entering Eq.(V.4) involve the derivatives and therefore do not vanish when G a = 0, so the 
variables £ t (£)are at variances 0(h 2 ) with the variables £ on the constraint submanifold. Respectively, MO and 
/(£) are at a variance 0(h 2 ) on the constraint submanifold also. 

The lowest order quantum corrections in Eq.(V.5) can originate from the Moyal bracket inside of the largest group 
of terms inside of the round brackets. In such a case, the symmetrized o-product can be replaced with the desired 
accuracy by the pointwise product. On the constraint submanifold, those terms do not contribute to the result and, 
respectively, the Moyal bracket in Eq.(V.5) can be replaced with the Poisson bracket. The second possible source 
of the quantum corrections is the operation o itself, as defined by Eq.(IV.15). To the lowest order in the Planck's 
constant, the fourth-order differential operator P 2 appears. The multiple product (...) o G ai ° G a2 --- ° G ak can be 
calculated as (...(((...) o G ai ) o G a2 )...) o G ak . The zeroth and second powers of H show up in the first three terms of 
the fc-series only. One gets, on the constraint submanifold, 

MO - {f(0,G^}V 2 G ai -^-^{{f(0,G^},G^}G ai V 2 G a2 

-^{{{f(0,G^},G^},G^}G ai G a2 r 2 G as +O(^). (V.7) 

The operator V acts on the left and on the right as prescribed by Eq.(IV.15). The terms G a involving no V action 
vanish on the constraint submanifold and are therefore omitted. The projected functions differ in general from the 
original ones. The constraint functions are, however, an exception, i.e., (G a )t = on the constraint submanifold and, 
furthermore, in the unconstrained phase space. 

The equivalence relations between the operators Eq.(III.7) lead, under the Weyl's association rule, to the equivalence 
relations between functions in the phase space: 

f(0 ~ 9(0 - MO = 9t(0 (V.8) 
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The physical observables are mapped through the quantization procedure onto the equivalence classes of functions. 
Since /(£) ~ ./t(£) and /(£) ^ f t (£,) on the constraint submanifold, the symbols ~ and « acquire distinct meaning. 
The average value of a function /(£) is defined as 

</>=/ (|^(0*/*(0*wi(fl ( v - 9 ) 

where is the symbol of the projection operator *P and W(£) is the Wigner function, i.e., the symbol of the density 
matrix r. In the classical limit 

P(ft 2m 

^o(2^ = n^(0)- (V.10) 

If the quantization is done in the standard canonical coordinate system, / '(2nh) m turns to the product of the 
delta functions to all orders in h. 

The evolution equation which is analogue of Eqs.(II.16) and (III. 16) takes the form 

|/(0 = /(Oa^(0- (v.ii) 

The quantum deformation of the Dirac bracket represents a skew-symmetric multiplication operation on the set of 
equivalence classes of functions in the phase space of quantum systems. 

Like in the classical case, the sum of two equivalence classes Ef and E g is defined by Ef+E g = Ef +g , the associative 
product EfEg is defined for quantum systems using the star-product as Ef titgt , and the skew-symmetric Dirac bracket 
is defined as {Ef, E g }o — Ef tAgt . These operations satisfy the Leibniz' law and the Jacoby identity. The associative 
product and the Dirac bracket do not depend on the choice of representatives of the equivalence classes. Since 
(G a ) t = 0, the deformed Dirac bracket vanishes for any two equivalence classes the one of which contains a constraint 
function. 

By quantization, the Dirac bracket is associated to the commutator of two projected operators. The Moyal bracket 
for projected functions on the constraint submanifold appears as a Lie bracket in the vector space O of equivalence 
classes Ef of functions in the phase space of quantum systems. The associative star-product and the skew- symmetric 
Moyal bracket gifts O with a Poisson algebra structure. 

The mapping Ef <-» f t induces further a Poisson algebra structure on the vector space of projected functions. The 
sum Ef + E g converts to ft + gt, the associative product EfE g converts to ft*gt, whereas the Dirac bracket {Ef, E g }i) 
becomes the Moyal bracket f t A g t . It is clear that these operations satisfy the Leibniz' law and the Jacoby identity 
and, furthermore, keep the vector space of projected functions closed, as (ft + gt)t — ft + gt, {ft * gt)t — ft* gt, and 
(ft Ah t ) t = f t Ah t . 

The skew-gradient projection of Eq.(V.ll) gives an evolution equation in the form involving projected functions 
only. Eq.(V.ll) is, however, more convenient for applications, since it is valid for any representative /(£) of an 
equivalence class and does not presuppose restrictions for /(£). 

In terms of the projected functions, the classical limit appears according to equations 

puft-kgt = fs9s, (V.12) 
h—>o 

lim/ t A fft = {f s ,9s} = {f,9}D, (V.13) 
h—>o 

where the use is made of Eqs. (11.15), (IV. 17) and (V.6). It is worthwhile to notice that Eq.(II.15) is valid on the 
constraint submanifold, so the operation f t A g t defined for all functions in the unconstrained phase space converts 
to the Dirac bracket on the constraint submanifold only. Away from the constraint submanifold, functions and, 
accordingly, the Dirac bracket do not make any physical sense. The equivalence classes and the projected functions 
are the only objects associated to physical observables. 
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In terms of the equivalence classes, the classical limit appears as 

VmEfEg = £ f £ g , (V.14) 
\im{E f ,E g } D = {£ f ,£ g } D . (V.15) 

The equivalence classes Ef of the quantum systems represent the quantum deformation of the equivalence classes £f 
of the corresponding classical systems. At h — ► 0, the star-product quantization recovers the classical constrained 
dynamics. 

To illustrate calculation of the quantum corrections to the Dirac bracket, in Appendix C we derive the analytical 
expression for the ti 2 correction. 



VI. QUANTUM SPHERICAL PENDULUM IN THE PHASE SPACE 

As an application we consider, in the phase space, the evolution of the Wigner function of a mathematical pendulum 
on an S n ~ 1 sphere of a unit radius in an n-dimensional Euclidean space with the coordinates <j) a . 

At the classical level, there exists, within the generalized Hamiltonian framework, two constraint functions 

g 1 =\n.4>, g 2 =n a <f) a (VIA) 

where = (0 a Q ) 1 / 2 ■ The first constraint Q 1 = implies that the particle stays on the sphere 0=1, the second 
constraint Q 2 = suggests that the radial component of the momenta vanishes. The constraint functions constitute 
a canonical pair (II. 2). One can check that the projected canonical variables £ s (£) = (0" , nf), with 

C = 070, (vi.2) 

7r s " = 07T" - Q 07T/0, (VI.3) 

and the Hamiltonian projected onto the constraint submanifold (see [24], Sect. 6), 

H s = l -{4>H^ -4> a 0VV 3 , (vi.4) 

are identically in involution with Q a , which is on a par with Eqs.(II.4) and (11.10). 

Since Q 2 is a second order polynomial with respect to the canonical variables, the infinite power series in the Poisson 
operator in Eq.(IV.16) is truncated at 0(h a ). Indeed, all the higher order terms vanish, so that Q 1 AG 2 = {G 1 1 Q 2 } = 1. 
Assuming G a = G a , Eqs.(V.2) hold. 

Similarly, Eqs.(V.3) hold for = £«(£) : G 1 and 0" depend on 0", so 0" AG 1 = 0. 7r" is a first degree polynomial 
with respect to the canonical momenta. In the expression 7r" A Q 1 , the operators J> 2k+1 act to the right on the 
coordinates only, and to the left on the momenta. The power series terms entering Eq.(IV.16) vanish starting with 
0(h 2 ), so irf A G 1 — {nf ,G 1 } = holds in virtue of Eqs.(II.4). G 2 is a second degree polynomial with respect to the 
canonical variables. The infinite power series in Eq.(IV.16) is truncated at 0(h°) again. Due to Eqs.(II.4), we have 
G 2 A — {G 2 , 0?} = and G 2 A < = {G 2 , <} = 0. 

Finally, one has to check that G a AH S = 0. The first equation, for a = 1, is valid since G 1 depends on 0", while 
H s depends on 7r" quadratically. The infinite power series in Eq.(IV.16) is truncated at O(h ), so the Moyal bracket 
can be replaced by the Poisson bracket. Since {G 1 1 7~L S } — 0, equation G 1 A TL S = holds. The second equation, for 
a = 2, is valid since G 2 is a second degree polynomial. Consequently, G 2 A H s = {G 2 ,H S } = 0. 

The projected canonical variables £ s (£), the constraint functions G a , and the classical Hamiltonian 7i s coincide with 
the projected canonical variables the constraint functions G a , and the Hamiltonian H t , respectively. 
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TABLE I: Skew-symmetric multiplication operations are listed for functions (second column) and equivalence classes of functions 
(third column) in the phase space of classical systems (first row) and in the phase space of quantum systems (second row) . 



Systems: 


unconstrained 


constrained 


classical 


{/,<?} 


U,9}d 


quantum 


/A 5 


ftf\gt 



The evolution equation for the Wigner function has the form (V.ll) where the sign of the right side should be 
changed. The power series expansion of the Moyal bracket over the Poisson operator is truncated at 0(h 2 ), since the 
Hamiltonian (VI. 4) is a fourth degree polynomial of the canonical variables, so we obtain 

8 h 2 ffiW 8 3 W 

-W = - { W M + - - 5^) - ^^(2^ - ^ ~ S»*»)). (VI.5) 

The first term in the right side is of the classical origin, while the second term represents a quantum correction to the 
classical Liouville equation and there are no other quantum corrections. Given IV (£,0) in the unconstrained phase 
space, W(£,t) can be found by solving the partial differential equation (VI.5). 

VII. CONCLUSION 

Real functions in the unconstrained phase space of second-class constraints systems split into equivalence classes 
corresponding to different physical observables. The quantum observables are described by equivalence classes of 
operators acting on the Hilbert space of states. The Weyl's association rule extends the equivalence relations to 
functions in the phase space of quantum systems. 

The Dirac bracket can be calculated as the Poisson bracket between functions projected onto the constraint sub- 
manifold using the phase flows generated by the constraint functions. The quantum deformation of the Dirac bracket 
is the Moyal bracket calculated for functions projected onto the constraint submanifold of the phase space. 

The skew-symmetric multiplication operations are synthesized in Table I. 

The operation f t A g t designates the quantum deformation of the Dirac bracket consistent with the canonical quan- 
tization of constraint systems. 

As an application of the general formalism, in Sect. VI we derived with the use of the operation f t A g t an 
evolution equation for the Wigner function of an n — 1-dimensional spherical pendulum, which represents a mechanical 
counterpart of the 0(n) non-linear sigma model. 

Our final conclusion states that quantum dynamics for constraint systems can be formulated in the unconstrained 
phase space provided symplcctic basis for constraint functions exists globally. 
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APPENDIX A: PROPERTIES OF THE SKEW-GRADIENT PROJECTION 

Let us apply Eq.(II.7) to a function /(£) expanded in a power series in £: 



3 



yl d p /(o) yl v — - i fif' 1 o ai \ g a2 \ G a ^\ 

p=0 r S ■> fc=0 fci + ...+fc-fc 1 P ~* 

fcl 

...{...{{e P ,a afc -^ +i },a afc -^ +2 },...g ofe K 1 ^ 2 -^ 



p=0 ^ s s fei>0,...,fc„>0 p ^ 

fcl 



...{...{{^,g a *-«^},g ak -«^},...G ak }Ga 1 Ga 2 -G, 



f'k 



= /(6(0)- (A.i) 

In this expression, the Taylor expansion is made first around £ = 0. To reach the step 3 of Eq.(A.l), the multibinomial 
formula for calculating the Poisson brackets of a product of p functions is used. Going from 3 to 4, the restriction 
k\ + ... + k p — k is removed. The summation over hi, ...,k p becomes thereby independent. To achieve the step 5 of 
Eq.(A.l), the summations over k i , . . . , kp are performed, which turns £ into £ s (£). The final result is given by Eq.(II.8). 
The classical counterpart of Eq.(III.6) follows straightforwardly from Eq.(II.8): 

(fi-f P U0 = hs(0-f PS (0- (A.2) 

In general, operators do not obey this property. However, under the circumstances specified by Eq.(III.6), this is still 
valid. Eq.(A.2) implies that 2 P products /i (£)•••/?(£) where /»(£) = /»(£) or = fi s (0 for 1 < i < p belong to the 
same equivalence class. Eq.(A.2) shows that class of the projected functions is closed under the pointwise product. 
Note other useful relation 

{f,9s}s = {fs,g s }. (A.3) 

To check it, we apply Eq.(II.7) for {/, g s } and use the Jacoby identity. Since g s is identically in involution with G a , 
one gets 

OC 

{.f,9s}s = Y,^ii-iif^ a ^^ a2 }^-^ ak }^s}Ga 1 Ga 2 ...Ga k 

k=0 
oc 1 

= E^{{-{{/^ oi }^ oa }.-^ *}^-i^-^,fl.} 

= {fs,9s}- 

In general, however, {/, g} s ^ {fs,g s }- To see this, one can set / = G a and g = G b ■ 
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The classical analogue of Eq.(III.17) is given by 

({f,9}D) s = {fs,9s} D . (A.4) 

This equation is a consequence of Eqs.(II.15) and, in particular, of the fact that the Leibniz' law applies to the Dirac 
bracket, 

{{/, 9}d, G a } = {{/, G a }, 9 }d + {/, {9, Q a }} D . (A.5) 

Using (III. 14), we get 

00 

({f,9} D )s = Y,^-mf^}D,Q ai },Q a2 },-Q ak }Ga i ga 2 -Ga k 

k=0 

= Etj E ^{{■■A{f^ ai }^ a2 },---^}A--\{9,Q ak ^ 1 },Q ak ^ 

k\ k 2 

Taking into account {,9 a }D — 0, one can place the last k constraint functions inside of the Dirac bracket. The 
summation over k can be removed and we obtain 

E ^{{■■\{f^ ai }^ a2 }i---G ak ^G*^--^ 

= {fs,9s}D- 

On the constraint submanifold, according to Eq.(II.15), the subscript can be omitted. In the unconstrained phase 
space, according to Eqs.(II.14) and the fact that f s and g s are identically in involution with the constraint functions, 
the subscript can be omitted also. Eq.(A.4) shows that class of the projected functions is closed under the Dirac 
bracket. 



APPENDIX B: PROPERTIES OF THE OPERATOR FUNCTION <8(j?) 

In order to check Eqs.(IV.3) - (IV. 8), it is useful to derive first the similar properties for 3S(r/): 

<8(t7)+ = (B.l) 
Tr [»(»?)] = {2Trh) n S 2n (r ] ), (B.2) 
®(0) - 1, (B.3) 

^-% {ri )Tr[%{-m = f , (B.4) 

Tr[%(ri)ib(-ri')] = (2Trh) n S 2n (r 1 - 7/), (B.5) 

»(»?)»(-»/) = K{r,-r,')eM-^lkilll kl )- (B.6) 

Eqs.(B.l) and (B.3) are conspicuous. Eq.(B.6) can be obtained using the identity e A+B = e A e B e~^ A ' B ^ that holds 
for operators A and B whose commutator is a c-number. Eqs.(B.2) and (B.5) can be derived by taking into account 
the explicit form of the matrix elements of 03 (77), 

n 

< >= S n (4>J -<% + Vn+1 ) exp(^- Va(<Pi + 4>2 )), (B.7) 

which can be obtained with the help of equation 

n 

exp(^ Vn+ a f l+a )\f< >= W - Vn+j > ■ 

a=l 
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Using Eq.(B.7) one gets 

r rj 2n n 

J J2^T < 0l|<8(r?)l<fc >K ^^(-^l^ >= S n (<f>? - <t>t)5 n {4>2 - 0?)- (B.8) 

Since Eq.(B.4) is valid for any operator, Eqs.(B.4) and (B.8) are equivalent. 

Applying the Fourier transform to 25(77), entering Eqs.(B.l) - (B.6), one gets Eqs.(IV.3) - (IV. 8). 
Let us multiply Eq.(IV.5) by fg, take the trace, and use Eq.(IV.12). As a consequence, we obtain 

™ = j$fcm°m. (b.9) 

The symmetrized star-product can be replaced with the pointwise product. Indeed, 

The quantum corrections 0(h p ) which could make the variance disappear, since they represent full derivatives and 
contribute to the surface integral only. By multiplying Eq.(IV.6) with g and taking the trace, one arrives at the same 
conclusion. 



APPENDIX C: LOWEST ORDER QUANTUM CORRECTION TO THE DIRAC BRACKET 

The series expansion over the Planck's constant is straightforward, so we restrict ourselves with the lowest order 
correction 0( ft 2 ). 

The quantum deformation of the Dirac bracket is associated to the operation 

St(OA/t(0 = (j(OA/.(0)t- (ci) 

Using Eq.(V.7), one can write 

flt(OA/t(0 = ff(0A/ t (0-f{{s(f)./.(0}.G B1 }P 2 G Bl 
~ ^{{{9(0Js(m,G^},G^}G ai V 2 G a2 

~^{{{{9(0Js(m,G^},G^},G^}G ai G a2 V 2 G a3 + O(h 4 ). (C.2) 

In terms 0{ ft 2 ), the Moyal bracket is replaced with the Poisson bracket and /t(£) is replaced with /«(£). The 
arguments similar to that used for the series expansion (V.7) allow to truncate the series expansion over G a , entering 



9t(0^ft(0 = ff(OA/t(0 

ft 2 .. ......... 1 

2! 

1 



—5-{{9(0, G ai }, f(0 + {/(£), G a2 }G Q2 + -{{/(£), G° 2 }, G a *}G a2 G a3 



3! {{{/(0. G a2 }, G° 3 }G a4 }G a2 G a3 G a4 }7> 2 G ai 

G" 1 }. G° 2 }, f(0 + {f(0,G a3 }G a3 + G a3 }, G a * }G Q3 G a4 }G ai V 2 G a2 

- G ai }, G" 2 }, G" 3 }, /(C) + {/(?), G^JGaJG,,,^? 2 ^ + 0(S 4 ). (C.3) 
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Let us turn to the first term of Eq.(C.3). We have 

ff(OA/*(0 = (s(Oa/(0). 



+ A G ai ) A G" 2 )... A G '- 1 ) A (<?(£) A G a <)) A G a ' +1 )... A G afc ) 



fc! 

fe=i ;=i 



oG ai O G Q2 ... O G afc + 
oo k 

+ E Ti £(-((/(0 A G ai ) A G" 2 )... A G a ") 



k\ 

fe=i ;=i 



3G ai O G a2 ... O Ga,^ O (#(£) A G a ,) O G a , + 1 ... O G ak . 



(C.4) 



The first term can be expanded using Eq.(V.7). In the second term, the Moyal bracket can be replaced with the 
Poisson bracket. The symmetrized star-product o is treated as in Eq.(V.7). The series over k to order 0(h 2 ) is 
truncated at k — 3 again. The third term is truncated at k = 4. The Dirac bracket originates, to the zeroth order in 
H, from the first term with k = and the third term with fc = 1 in Eq.(C.4). The result takes the form 

ff(OA/ t (0 = {$(0, f(0}n 
h 2 



n 2 l 



-^{{.f(0,9(0},G^}V 2 G ai - _-{{{/(^),^)},G ai },G Q2 }G ai ^G Q 



8 2! 

{{{{/(0. 9(0}, G ai }, G" 2 }, G a3 }G ai G a2 P 2 G a3 



-(f(^G^){g(0,G ai } _{/(e),G ai }(5(C)^G ai ) - -{/(£), G Ql }7> 2 {<?(0, G a J 



~~8~3! 

_ft 2 
24 

-y{/(0,{ff(0,G ai }}^ 2 G ai 
-f^{{/(0,{5(0,G ai }},G-}G ai P 2 G a 

8 2 
J?l 

8 3 
J? I 

8 3 
1 

8 3 

~ ^({{f(0,G^},G^}r 2 G ai ){g(0,G a2 } 



{{/(0,G ai },{ 5 (e),G a2 }}G ai P 2 G Q 



-5-^t{{{/(0. G 01 }}. G a2 }, G a3 }G ai G a2 'P 2 G a 



{{{fit), G ai }, {9(0, G° 2 }} , G a3 }G ai G a2 T' 2 Ga 
{{{/(£), G ai }, G a2 }, {«?(£), G a3 }}G ai G a2 P 2 G Q 



" 8 2 



ft 2 1 
8 

ft 2 1 



{{/(C),G ai },G a2 }G ai P 2 { 5 (C),G a2 } 



-y^{{/(C),G ai },G a2 }{ 5 (e),G ai }P 2 G a2 



({{{f(0,G ai }, G a2 }, G a3 }G ai P 2 G O2 ){ 5 (0, G a3 } 



-y^{{{/(0,G ai },G a2 },G a3 }G ai G a2 7' 2 {.9(0,G a3 } 

ft 2 2 
" 8 3 

~ G ai }, G a2 }, G a3 }, G°*}{ 5 (0, G ai }G a2 G a3 7 :,2 G a4 + 0(ft 4 ) 



{{{/(£), G ai }, G a2 }, G°n{5(0, G ai }G a2 P 2 G Q 



(C.5) 



In the first line, we recover the Dirac bracket (11.14). The second term is the quantum correction to the Poisson 
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bracket. The third and fourth lines originate from the expansion of the first term in Eq.(C.4) as prescribed by 
Eq.(V.7). The fifth line comes from the expansion of the k = 1 component of the third term in Eq.(C.4). The lines 6 
- 11 appear from the expansion of the second term in Eq.(C.4). The rest appears from the components k > 1 of the 
third term in Eq.(C.4). The round brackets unify terms within the action of the operation V 2 . If the round brackets 
are suppressed (in the most cases), V 2 acts on all terms. 

Combining Eqs.(C.3) and (C.5), we get the Dirac bracket on the constraint submanifold to the order 0(h 2 ). 
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